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HIGHLIGHTS 


►  An  equivalent  circuit  is  used  to  describe  the  characteristics  of  LiFeP04  batteries. 

►  Three  different  model-based  algorithms  are  designed  to  estimate  the  battery  SOC. 

►  The  estimation  approaches  are  verified  with  two  typical  driving  profiles. 

►  The  system  robustness  and  the  convergence  behavior  of  SOC  estimators  are  compared. 
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One  of  the  most  important  aspects  in  battery  management  systems  (BMS)  in  electric  vehicles  is  the  state 
of  charge  (SOC)  estimation.  SOC  needs  to  be  accurately  determined  for  safety  and  performance  reasons 
but  cannot  be  measured  directly  due  to  the  flatness  and  hysteresis  of  the  open  circuit  voltage  (OCV)  curve 
of  Lithium-ion  chemistries  as  LiFeP04.  The  classical  approach  of  current  integration  (Coulomb  counting) 
cannot  solve  the  problems  of  accumulative  error  and  inaccurate  initial  values,  thus  advanced  estimation 
algorithms  are  applied  to  determine  the  sate  of  charge.  In  this  work,  three  model-based  state  observer 
designs  including  Luenberger  observer.  Extended  Kalman  Filter  (EKF)  and  Sigma  Point  Kalman  Filter 
(SPKF)  are  carried  out  and  studied.  These  estimation  approaches  are  verified  using  measurement  data 
acquired  from  commercial  LiFeP04  cells.  In  addition,  computational  tests  analyze  the  systems  perfor¬ 
mances  in  terms  of  tracking  accuracy,  estimation  robustness  against  temperature  uncertainty,  sensor 
drift,  and  convergence  behavior  with  an  initial  SOC  offset. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Recently,  with  more  emphasis  on  environment-friendly  vehi¬ 
cles,  the  use  of  electric  powered  transportation  is  increasing.  For 
automotive  industry  the  lithium-ion  energy  storage  system  is 
currently  considered  as  the  best  alternative  to  meet  high  specific 
energy  and  specific  power  demands  of  electric  vehicles  (EVs).  Based 
on  current  research  and  development,  LiFeP04  (LFP)  as  active 
material  of  the  positive  electrode  is  more  stable  and  safer  as  other 
cell  chemistries  under  abusive  operating  conditions  [1],  Together 
with  its  low  cost  and  long  lifespan,  LFP  is  a  favorable  choice  as 
cathode  material  in  EV  applications. 

The  battery  pack  of  an  EV  is  controlled  by  a  battery 
management  system  (BMS).  In  complex  and  dynamic  vehicle 
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operation  conditions,  a  good  BMS  is  required  to  provide  accurate 
knowledge  of  the  battery’s  state  of  charge  (SOC).  This  infor¬ 
mation  reflects  the  remaining  capacity  that  can  be  drawn  from 
the  battery  pack  and  is  used  to  ensure  an  optimum  control  of 
charging  and  discharging  processes  [2],  Flowever,  several  factors 
for  SOC  determination  for  LFP-batteries  such  as  hysteresis 
phenomena,  flat  characteristic  of  open-circuit-voltage  (OCV)  over 
SOC,  and  limited  voltage  measurement  accuracy  can  result  in  SOC 
estimation  errors. 

A  widely  used  method  to  estimate  SOC  is  the  Coulomb 
counting  approach  [3,4],  which  samples  the  battery  current 
and  computes  the  accumulated  charge  and  discharge  to  estimate 
SOC.  This  approach  is  simple  and  inexpensive  to  implement. 
However,  it  faces  drawbacks  such  as  accumulated  measurement 
error  and  inaccurate  initial  SOC  values.  Other  methods 
described  in  [5],  such  as  discharge  test  or  internal  resistant 
characteristic,  are  usually  considered  as  laboratory  methods. 
Additionally,  the  artificial  neural  networks,  e.g.  [6]  can  simulate 
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the  complicated  battery  dynamic.  However,  its  performance  is 
sensitive  to  training  data.  Therefore,  it  is  not  appropriate  for  the 
battery  online  application. 

Further  developed  strategies  for  battery  SOC  estimation  are 
model-based  techniques.  State  estimation  with  Luenberger 
observer  has  already  been  analyzed  [7],  Later  on,  a  Kalman  Filter 
(KF)  approach  with  its  knowledge  about  statistical  characteristics  of 
process  and  measurement  noise  is  intensively  studied  [8-10]. 
Extended  Kalman  Filter  (EKF)  and  Sigma  Point  Kalman  Filter  (SPKF) 
are  two  improvements  to  the  classical  KF  which  are  investigated  for 
nonlinear  systems.  EKF  uses  the  first  order  Taylor  series  to  trans¬ 
form  the  nonlinear  system  into  a  linear  system  [11]  while  the 
second  approach  uses  the  unscented  transformation  to  approxi¬ 
mate  the  probability  density  function  of  the  nonlinear  system  with 
series  deterministic  points  [12], 

This  work  presents  a  survey  and  a  comparative  study  of  the 
model-based  SOC  estimation  approaches  above  with  focus  on 
estimation  robustness.  At  first,  the  considered  estimation  algo¬ 
rithms  are  separately  studied.  Based  on  the  theoretical  and  prac¬ 
tical  analysis,  this  work  determines  the  optimal  choice  to  estimate 
dynamic  SOC  behavior  for  LFP  cells. 

2.  Description  of  cell  model 

A  battery  equivalent  circuit  used  for  model-based  state  esti¬ 
mation  design  is  shown  in  Fig.  1. 

The  steady  state  cell  behavior  is  characterized  by  the  OCV 
voltage.  It  refers  to  the  potential  difference  between  two  electrodes 
in  an  open  circuit.  At  a  certain  temperature,  the  relationship 
between  SOC  and  OCV  is  obtained  through  measurement  and  can 
be  established  for  the  state  estimation  using  a  nonlinear  function  or 
lookup  tables  [11,13],  The  dynamic  cell  behavior  is  described  by  an 
ohmic  resistance  Ro  with  two  RC  elements  in  series  with  a  short  and 
a  long  time  constant.  The  parameters  can  be  identified  using 
nonlinear  least-squares  methods,  which  minimize  the  error 
between  the  model  output  and  experimental  data  until  a  pre¬ 
defined  error  bound  is  reached  [15,16].  Based  on  the  Coulomb 
counting  law  and  the  linear  differential  equation  of  RC  circuits,  the 
discrete-time  dynamic  equations  of  the  cell  model  and  the  cell 
voltage  can  be  expressed  by 

SOC(k)  =  SOC(k-  1)  +/(k)-?7-At/CN,  (1) 

Vt(k)  =  (l-At/^COJ.V^k-lJ  +  ^kJ-At/Ct,  (2) 

and 

V2(k)  =  (1  -  \t/(R2C2))-V2(k-  1)  +/(/<) -At/C2,  (3) 


resistances  of  RC  circuits  Ri  and  R2  as  well  as  capacitances  of  RC  circuits  C]  and  C2. 


where  k  is  the  discrete-time  index,  Vi  and  V2  the  capacitor  voltages, 
t]  the  columbic  efficiency,  and  Cn  the  discharge  capacity  of  the 
battery.  Additionally,  charging  current  is  defined  as  a  positive  value. 
Hence,  the  cell  voltage  is  calculated  as 

Vcell(k)  =  OCV(SOC  (fc))  +  V,  (k)  +  V2(k)  +  R0I(k).  (4) 

By  referring  to  Fig.  1  and  considering  the  current  I  as  the  model 
input,  the  state  vector  of  the  cell  model  xk  can  be  chosen  as 

X  =  [SOC,Vi,V2]t.  (5) 

In  order  to  facilitate  the  application  of  model-based  estimation 
algorithms,  the  system  equations  from  (1)  to  (4)  need  to  be  trans¬ 
lated  into  a  discrete-time  state  space  form 

xfc  =  Axk  i  +  B  uk,  x(t0)  =  Xq,  (6) 

Zic  =  Vcell  =  ^(Xfc^.Ufc),  (7) 

where  A  is  the  dynamic  matrix,  B  the  input  matrix,  g(  )  the 
nonlinear  process  output  model,  u  the  model  input,  and  xo  the 
initial  state  at  time  t  =  to. 

3.  Proposed  methodology 

Three  closed  loop  state  observers  are  investigated  in  this  work. 
The  general  structure  of  battery  SOC  estimation  is  shown  in  Fig.  2. 

The  principle  of  SOC  observation  is  to  generate  the  cell  voltage 
from  the  state  estimator  and  to  compare  it  with  the  measurement 
value.  By  feeding  back  this  difference  into  the  model,  the  predicted 
battery  SOC  can  be  corrected  with  a  so-called  feedback  gain  matrix. 
This  matrix  can  be  determined  by  the  following  methods: 

3.1.  Luenberger  observer 

The  Luenberger  observer  is  widespread  for  linear  deterministic 
dynamic  systems  in  the  industry  because  of  its  simple  imple¬ 
mentation  and  tuning  of  the  observer  gain.  Based  on  the  system 
knowledge,  the  Luenberger  observer  is  usually  designed  with 
a  constant  feedback  gain.  Therefore,  in  spite  of  a  nonlinear  output 
function  (4),  the  state  correction  is  proportional  to  the  voltage 
residual  ek 


Fig.  2.  The  structure  of  model  based  battery  SOC  estimation  with  cell  temperature  T, 
terminal  voltage  Vcen,  input  u,  real  states  x,  and  estimated  states  x. 
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ek  =  zk  ~  Zfc, 

(8) 

and 

xk  =  Axm  +  Buk  +  kek, 

(9) 

with  the  Luenberger  gain  vector  k. 

Note  that  the  OCV  of  LFP  cells  is  very  flat  over  a  wide  SOC  range. 
The  residual  ek  calculated  between  measured  and  modeled  cell 
voltage  may  not  correspond  to  the  actual  SOC  variation  due  to  the 
measurement  noise.  Therefore,  two  Kalman  based  filter  algorithms 
are  further  used  for  this  nonlinear  estimation  problem. 

3.2.  Extended  Kalman  filter 

The  Kalman  filtering  is  an  optimal  recursive  solution  to  the 
linear  Gaussian  problem.  A  Kalman  filter  estimates  the  state  of 
a  dynamic  system  with  a  linear  process  model  and  measurement 
model  [17],  To  nonlinear  problems,  the  most  well  known  inference 
of  the  Kalman  Filter  framework  is  the  extended  Kalman  filter  (EKF). 
Compared  to  a  linear  state  space  description  in  (6),  the  general 
framework  of  a  discrete-time  nonlinear  system  with  noise  input  is 
given  by 

Xfc  =  f(xk_i,  ufe} +  wfc,  x(t0)  =  Xo,  (10) 

2k  =g(Xfc_l,Ufc)  +  Vk,  (11) 

where  f(  )  and  g(  )  are  the  nonlinear  process  model  and  the 
measurement  model.  Vector  wk  represents  the  system  noise  and  vk 
the  measurement  noise.  Both  wk  and  vk  are  assumed  to  be  uncor¬ 
related,  zero-mean  Gaussian  white  sequences  and  their  covariances 
are  known  as  Qk  and  rfe  respectively  [17],  The  EKF  algorithm  could 
be  mathematically  expressed  as 

•  Prediction  update 

Xfc  =  f(xfc_i,ufc,Xo),  (12) 


Notice  that  the  key  idea  underlying  the  EKF  approach  is  the 
linearization.  It  approximates  the  nonlinear  function  utilizing  the 
first  order  Taylor  expansion  by  the  partial  derivative.  A  detailed 
explanation  for  EKF  is  given  in  [18],  Accordingly,  one  important 
limitation  of  the  EKF  is  its  local  linear  expansions.  For  the  battery 
SOC  estimation,  because  of  the  stepwise  flatness  and  the  hysteresis 
effect  of  the  OCV  and  SOC  curve  the  measurement  model  is 
inherently  nonlinear.  When  the  degrees  of  system  uncertainty  or 
nonlinearity  are  significant,  these  approximations  of  the  EKF  can 
result  in  large  estimation  errors  [12]. 

3.3.  Sigma  point  Kalman  filter 

The  Sigma  Point  Kalman  Filter  (SPKF)  or  unscented  Kalman 
Filter  (UKF)  addresses  the  approximation  issues  of  the  EKF.  Fig.  3 
illustrates  the  linearization  applied  by  the  sampling  methods,  EKF 
and  SPKF.  Instead  of  local  linearization,  the  SPKF  captures  the 
statistical  distribution  of  the  nonlinear  system  using  deterministic 
sampling  (so-called  sigma  points). 

In  the  general  case,  these  sigma  points  are  located  at  the  mean 
and  symmetrically  along  the  main  axes  of  the  covariance  [18],  The 
unscented  Kalman  Filter  can  be  summarized  as  follows 

•  Calculating  sigma  points  with  mean  x&_i  and  covariance  P/( 
X[0]  =  **-!,  (18) 

XW  =x*-1  +  (V(n  +  A)P*)  i  =  1,2,  ...,n,  (19) 

Xw  =Xfc-i-(vW*)pk)  ,  i  =  n  +  l,n  +  2,...,2n,  (20) 

where  n  is  the  dimension  of  the  state  vector.  The  number  of  sigma 
points  x  is  calculated  as  2n  +  1.  The  parameter  A  is  expressed  as 

A  =  a2(n  +  k)  -n,  (21) 


Pfc  =  AkPfc  ,A[  +  wkQk  ,w[.  (13) 

and 

•  Measurement  update 

kk  =  PkCTk(ckPkCTk+vkr  kvTky\  (14) 

xk  =  xk  +kk(zk-g(xk,uk)),  (15) 

Pk  =  (l-kkCk)Pk,  (16) 


where  Pk  is  the  error  covariance  of  the  state  estimation,  xk  the 
priori  state  estimate  at  step  k  and  xk  the  posterior  state  estimate. 
Ak,  Bk  are  defined  as  the  Jacobian  matrices 


Ak 


9g| 

;tx|x='  ' 


(17) 


(a)  Actual  (b)  Linearized  (EKF)  (c)UT(SPKF) 


sigma  point  approach  [18], 
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where  a  (0  <  a  <  1)  and  k  (default  choice  is  k  =  0)  are  tuning 
parameters  that  determine  the  distribution  of  the  sigma  points 
spreaded  from  the  mean  [19]. 

•  Calculating  weights  associated  with  each  sigma  point 


W™  =  ^TI+(1-a2  +  /?)’  (23) 

w™  =  w^1  «  1  i  =  1.2.  ...,2n.  (24) 

2n  +  2A  v 


where 


ii 

£ 

(32) 

p,  ,  -*,)(?? -ik)+vkvl 

2  =  0 

(33) 

and 

p=  =  l:»p(*?-*‘)(2*"1-i>)r 

2  =  0 

(34) 

4.  Experimental  details 

where  wm  and  wc  are  sets  of  weighting  matrices.  The  parameter 
/?  can  be  used  to  incorporate  the  higher  order  knowledge  about  the 
distribution.  For  a  Gaussian  distribution,  =  2  is  the  optimal  choice 

(191. 

•  Prediction  update 

Update  each  sigma  point  through  the  nonlinear  dynamic 
equation  to  predict  the  system  state  as 

Xk  =  f  (Xk-i,Uk),  (25) 

(26) 

2  =  0 
and 

Pk  =  (C"  -  X*)  (x7]  -  *k)  +wkQ.k  iK  (27) 

i  =  0 

where  the  prior  state  estimation  xk  is  derived  from  the  weighted 
mean  of  the  sigma  points  vector  xk.  The  prior  covariance  matrix  Pfc 
is  estimated  as  the  weighted  sum  of  the  squared  distances  of  all  the 
sigma  points  from  xk . 

•  Measurement  update 

Update  each  sigma  point  through  the  nonlinear  measurement 
equation  to  predict  the  measurement  value  as 

Zk=  gfxkMk);  (28) 

%  =  f>£zE  (29) 

i=o 

After  the  prediction  step,  the  estimates  of  the  system  states  and 
covariance  are  defined  by 

xk  =  +kfc(zfc-ifc),  (30) 

and 


In  the  present  work,  an  LFP  pouch  cell  with  high  specific  energy, 
having  a  nominal  capacity  of  20  Ah  and  a  nominal  voltage  of  3.3  V, 
is  used  for  the  comparison  study  of  the  SOC  estimation  algorithms. 
The  SOC  verification  tests  are  carried  out  on  a  BaSyTec  test  bench 
with  0-10  V  voltage  range  and  currents  up  to  250  A.  The  accuracy 
of  both  voltage  and  current  measurement  is  within  0.05%  of  the  full 
scale  range,  and  the  resolution  within  0.003%.  During  the  test,  the 
cell  is  located  in  a  climatic  chamber  ESPEC  PU-3KP  (temperature 
range  -40  °C-100  °C,  accuracy:  ±0.5  °C)  and  the  cell  ambient 
temperature  is  constantly  held  at  30  °C. 

In  order  to  evaluate  the  SOC  estimation  approaches,  two  typical 
driving  profiles,  the  New  European  Driving  Cycle  (NEDC)  and  the 
ARTEMIS  cycle  [14],  are  applied.  The  current  sequences  are  derived 
from  these  driving  profiles  and  depicted  in  Fig.  4.  The  first  test  with 
one  NEDC  current  load  takes  approximately  20  min.  The  initial  SOC 
is  adjusted  to  80%  SOC.  In  the  second  test,  the  entire  operating 
range  of  the  cell  is  utilized.  The  ARTEMIS  load  cycle,  started  with 
90%  SOC,  is  repeated  until  the  cell  voltage  reaches  its  limitation 
(1  min  rest  between  two  cycles).  Data  points  are  collected  at  50  Hz 
during  the  tests. 

The  nonlinear  battery  model  in  Matlab/Simulink  incorporates 
the  measured  data  including  voltage,  current,  and  temperature 
from  the  test  bench  and  estimates  the  SOC  value  with  the  presented 
approaches.  The  SOC  performances  are  compared  with  the  current 
integration  which  serves  as  the  actual  SOC  behavior. 

5.  Results  and  discussion 

5.1.  NEDC  with  80%  start  SOC 

Figs.  5  and  6  demonstrate  the  results  of  the  proposed  estima¬ 
tion  approaches  with  the  NEDC,  when  the  SOC  is  correctly 
initialized  to  80%.  Tracking  for  SOC  behavior  is  shown  in  Fig.  5. 
Fig.  6  illustrates  the  cell  terminal  voltage  estimation.  Generally,  the 
estimated  value  is  consistent  with  the  measured  one.  Nevertheless, 
using  nonlinear  algorithms,  SOC  can  be  estimated  more  accurately 
(error  less  than  1%)  in  comparison  with  2%  error  with  the  Luen- 
berger  observer. 

Note  that  the  estimation  errors  start  at  zero  and  diverge  slowly 
away  from  the  true  behavior.  Although  the  voltage  estimation  error 
of  SPKF  is  approximately  zero  (±1  mV  in  Fig.  6),  the  uncertainty  of 
the  battery  model  could  result  in  the  slight  SOC  divergence. 

5.2.  ARTEMIS  driving  cycle  with  observable  start  SOC  and 
unobservable  start  SOC 


P  k 


=  Pk  -  kfcPzk[, 


(31) 


In  the  second  test  regime  the  SOC  algorithms  are  performed 
with  the  ARTEMIS  driving  cycle  including  a  more  dynamic  current 


Load  current  [A]  Load  current  [A] 
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Fig.  4.  Applied  current  sequences  of  NEDC  and  ARTEMIS  driving  profiles  utilized  for 
the  SOC  estimation. 


Time  [s] 


Time  [s] 

error  under  NEDC. 


load  for  the  entire  operation  range.  Fig.  7  depicts  the  estimation 
results  with  both  observable  and  unobserverable  start  SOC. 

The  achieved  SOC  performances  are  evaluated  by  the  maximum 
estimation  error  and  the  convergence  behavior  for  incorrect 
initialization  of  SOC.  When  the  start  SOC  is  correctly  initialized  to 
90%,  all  approaches  give  similar  estimation  errors  (maximum  3% 
error).  However,  an  observable  start  SOC  usually  cannot  be  guar¬ 
anteed  in  practice.  By  defining  a  20%  start  SOC  error,  the  superior 
performance  of  SPKF  can  be  clearly  seen.  The  SPKF  enables  the 
fastest  convergence  rate  and  stabilization  of  SOC  behavior  without 
remarkable  SOC  oscillations. 

5.3.  Temperature  uncertainty  and  sensor  drift 

Cell  temperature  is  one  of  the  most  significant  factors  impacting 
both  the  performance  and  the  useful  life  of  the  battery.  Because  of 
different  sensor  positions,  battery  aging  and  cooling  strategy  of 
BMS,  the  cell  temperature  could  not  be  directly  or  exactly 
measured.  In  order  to  investigate  the  estimation  robustness  of 
proposed  approaches,  a  +3  °C  temperature  offset  and  a  +4  mV 
voltage  measurement  offset  are  further  embedded  in  the  computer 
simulation. 

Simulations  are  carried  out  with  the  same  ARTEMIS  drive  profile 
starting  from  90%  initial  SOC  as  for  Fig.  7(a).  In  Fig.  8(a)  the  simu¬ 
lated  temperature  is  raised  with  3  °C,  which  reduces  the  model 
accuracy  and  changes  slightly  the  SOC-OCV  behavior.  Fig.  8(b) 
shows  the  result  with  a  voltage  measurement  error.  Very  close 
agreement  in  Fig.  8  is  observed,  as  expected.  It  can  be  concluded 
that  the  SOC  performance  with  Luenberger  observer  is  more 
influenced  by  the  measurement  accuracy.  The  estimation  error 
from  both  situations  is  enlarged  and  lasts  to  the  end  of  the  test.  In 
contrast,  behaviors  of  EKF  and  SPKF  are  nearly  consistent  with  the 
measurements  from  Fig.  7(c). 


Fig.  5.  Comparison  of  SOC  behavior  and  estimation 
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Fig.  7.  Comparison 


Estimation  error  under  temperature  offset 


and  +4  mV  sensor  drift  for  voltage  measurement  (b). 


Table  1 

Root  mean  square  error  (RMSE)  of  different  methods  in  SOC  estimation  for  ARTEMIS 
driving  cycle  from  90%  SOC. 

Methods  Observable  Unobservable  Temperature  Sensor 

start  SOC  start  SOC  uncertainty  drift 

Luenberger  1.05  7.3  1.8  2.7 

EKF  0.95  6.2  0.94  1.8 

SPKF  0.96  4.3  0.95  1.6 


To  summarize,  an  overview  of  root  mean  square  errors  (RMSE) 
for  each  estimation  method  in  terms  of  considered  situations  is 
presented  in  Table  1.  The  classical  Luenberger  observer  relies 
mostly  on  the  accuracy  of  the  battery  model,  as  it  can  hardly  take 
the  model  uncertainty  and  measurement  noise  into  account.  EKF  is 
comparatively  more  accurate  than  Luenberger  method.  The 
difference  between  EKF  and  SPKF  is  mainly  reflected  in  the 
convergence  rate  with  an  unobservable  start  SOC.  For  most  states, 
SPKF  provides  the  better  SOC  estimation  results  than  other 
approaches. 

6.  Conclusion 

In  this  work,  three  model-based  algorithms  including  standard 
Luenberger  observer,  EKF,  and  SPKF  are  applied  to  monitor  the  SOC 
of  an  LFP  battery.  First,  an  equivalent  electrical  model  is  represented 
to  describe  the  battery  electrical  behavior.  Then,  the  general 
formulations  of  estimation  approaches  have  been  studied.  Their 
performances  are  compared  using  the  measurement  data  acquired 
from  a  commercial  LFP  cell  on  a  test  bench.  Simulation  results  denote 
the  improvement  of  SPKF  for  the  battery  SOC  estimation  in  terms  of 
tracking  accuracy,  the  convergence  behavior,  and  the  estimation 
robustness.  In  addition,  the  nonlinear  approximation  using 
unscented  transformation  avoids  computing  the  Jacobian  matrix  of 
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EKF  (17).  Therefore,  SPKF  can  also  provide  a  better  numerical 
stability. 

In  future  work,  the  algorithms  will  be  implemented  on  a  BMS 
and  thereby  be  validated  under  limited  achievable  computation 
power  and  precision  of  floating  point  numbers.  In  a  further  devel¬ 
opment  stage  the  underlying  battery  model  will  be  extended  taking 
battery  aging  (capacity  fade  and  resistance  rise)  into  account. 
Therewith  the  performance  of  the  algorithms  over  the  lifetime  of 
the  battery  will  be  further  improved. 
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